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ABSTRACT 

Polarized emission from a quasar is produced by wavelength-independent electron 
scattering surrounding its accretion disc, and thus avoid the contamination from its 
host galaxy and reveal the intrinsic emission spectrum of the accretion disc. Ultra- 
violet (UV) emission from a quasar is normally free from the contamination from 
its host galaxy. Polarization fraction of the quasar's disc emission can therefore be 
determined by comparing total UV emission with polarized visible to near-infrared 
(NIR) emission; and the resulting continuum spectrum from UV to infrared can reveal 
the theoretically expected Balmer edge absorption feature. We fit the polarized spectra 
in visible and NIR bands together with the total UV spectra of two type-1 quasars 
(3C 95, 4C 09.72), to an extended geometrically thin and optically thick accretion disc 
model. In addition to the standard model, we include the Balmer edge absorption due 
to co-rotational neutral gas on a narrow annulus of the accretion disc. We find that the 
extended thin accretion disc model provides adequate description on the continuum 
spectra of the two quasars from UV to NIR wavelengths. A Monte-Carlo-Markov- 
Chain fitting to the continuum spectra is able to well constrain the true polarization 
fraction of the disk emission, which allows the Balmer edge feature to be completely 
revealed from polarized visible to UV continua. The Balmer edge feature is prominent 
in both quasars' spectra, and is significantly broadened due to the orbital motion of gas 
in the accretion disc. The broadening of the Balmer edge feature is therefore related 
to the quasar's inclination. This work proves the concept of determining quasar's 
inclination from the Balmer edge feature in their continuum spectra. 

Key words: galaxies: active - quasars: general - techniques: polarimetric - tech- 
niques: spectroscopic - methods: data analysis. 



1 INTRODUCTION 

Quasars are the most luminous objects in the Universe. The 
release of gravitational energy as material is accreted onto 
the supermassive black hole (SMBH) heats the accretion disc 
and generates multi-color thermal emission from ultraviolet 
(UV) to infrared (IR) wavelengths (e.g., Shakura & Sun- 
yaev 1973; Novikov & Thorne 1973; Lynden-Bell & Pringle 
1974; Laor & Netzer 1989; Laor et al., 1990; Hubeny et al. 
2000). However, the total emission of a quasar at visible 
and near-IR (NIR) wavelengths is mixed by emissions from 
both its accretion disc surrounding its central SMBH and 
its host galaxy, which complicates observations of intrinsic 
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emission of the quasars' accretion disc. For example, the op- 
tically thick and geometrically thin "standard" disc model 
predicts the emission spectrum in the optical to NIR should 
follow F v oc v a , where F u is the flux per frequency, v is 
the frequency and a = 1/3 (e.g., Shakura & Sunyaev 1973; 
Lynden-Bell & Pringle 1974) . However many quasar spectra 
in the optical to NIR have a < —0.2 (Neugebauer et al. 1987; 
Cristiani & Vio 1990; Francis et al. 1991; Zheng et al. 1997), 
consistent with significant host galaxy contaminations which 
are increasingly important at longer wavelengths. 

Polarized spectra of quasars at visible and NIR wave- 
lengths may reveal the emissions of their accretion discs, 
because polarization is expected to come from the electron 
scattering inside the Broad Line Region (BLR) and the elec- 
tron scattering is wavelength independent. Kishimoto et al. 
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(2008) reported polarized spectra of six quasars in the rest 
wavelength ranging from 0.2 - 2 pirn. These quasars have ~ 
1% of polarization in optical bands, less than 1% of polar- 
ization in NIR continuum, and little polarization in their 
emission lines. Indeed Kishimoto et al. (2008) found that 
the polarized NIR spectra manifest a = 1/3, expected from 
the standard disc model. We emphasize two important facts 
regarding quasars' polarization at visible and NIR wave- 
lengths: first, the alignment of optical polarization with the 
radio structure of many quasars suggests that the emission 
is polarized not in the disc atmosphere but in an equato- 
rial scattering region surrounding the disc (e.g., Stockman 
& Angel 1979; Smith et al. 2005). Indeed, intrinsic polariza- 
tion produced in the disc atmosphere is likely to be damped 
by Faraday rotation (e.g., Agol & Blaes 1996). Second, the 
observed polarization percentage in the optical to NIR bands 
may not be the polarization percentage of the intrinsic disc 
emission, or the polarization percentage of the equatorial 
scattering region, due to possible contaminations from the 
quasar's host galaxy. 

The standard thin disc model predicts a "bump" shape 
from UV to visible wavelengths, which is produced by the 
emission from the inner disc region and is determined by 
the combination of the mass of the central supermassive 
black hole (BH) and the inner disc boundary radius. The 
UV continuum is relatively free from the host galaxy's con- 
tamination, and thus should reflect the emergent spectra of 
the accretion disc and allow to probe the inner part of the 
accretion disc. Two quasars of the Kishimoto (2008) sample 
(3C 95, 4C 09.72) have been observed with the Hubble Space 
Telescope (HST) and have published UV spectra (Bahcall 
et al. 1993; Marziani et al. 1996; Evans & Koratkar 2004). 
Combining the UV continuum and the polarized spectra in 
visible and NIR wavelengths, therefore, will provide pow- 
erful spectral diagnosis of the nature of the two quasars' 
supermassive BHs and accretion discs. It has also been re- 
ported that the quasar continuum spectra in wavelengths 
shorter than 120 nm can be described by a simple power 
law with index o?euv = —1.76 ±0.12, which suggests that 
the extreme UV (EUV) continuum is probably due to the in- 
tergalactic medium photoionized by the integrated radiation 
from quasars (Telfer et al. 2002). Therefore the EUV con- 
tinuum may not reflect the intrinsic emission of a quasar's 
accretion disc. 

Quasars' continuum spectra from UV to NIR wave- 
lengths may be used to constrain the properties of their 
central supermassive BHs, in particular their BH spin, in 
a similar way to the method used for stellar mass BHs in 
X-ray binaries whose BH mass and disc inclination angle (i) 
are known very well for some of them (Zhang et al. 1997; 
Shafee et al. 2006). Czerny et al. (2011) used the standard 
thin accretion disc model to fit several broad-band photo- 
metric points of quasar SDSS J094533. 99+100950.1. Assum- 
ing a BH mass, this procedure put constraints on its spin. 
However, the broad-band photometric fluxes may be con- 
taminated by broad spectral lines, which may affect the es- 
timates of the BH spin. Davis & Laor (2011) also used the 
thin accretion disc model to fit the shape of optical contin- 
uum of a set of 80 PG quasars, and obtained their radiative 
efficiency and derived their BH spin. However, both meth- 
ods had to assume the BH mass and the inclination angle 
of the accretion disc a priori. Various methods can be used 



to estimate the BH mass in a quasar, but normally within 
a factor of 2-3, as will be discussed later in section 4.1. 

The inclination of the accretion disc for an individual 
quasar is normally not well constrained. According to the 
unified model of active galactic nuclei, the emission from 
a quasar's accretion disc is not blocked by its dusty torus, 
which implies a disc inclination of less than ~ 60° (e.g., An- 
tonucci 1993). Advanced torus model and spectroscopy of 
quasars in the mid-infrared wavelengths may provide more 
insights on the interaction of the disc emission and the 
torus and better constrain the disc inclination (e.g., Alonso- 
Herrero et al. 2011). Another method to determine the disc 
inclination is to describe not only the spectra but also the 
morphology of active galaxies self-consistently (Kacprzak et 
al. 2011). For radio loud objects, inference of radio core lumi- 
nosity may constrain the view angle of the radio source, i.e., 
the inclination of the accretion disc (e.g., Wills & Brotherton 
1995). However, the disc inclination has not been estimated 
from the disc emission directly. In principle, the luminosity 
of the accretion disc scales as cos(i), so the inclination may 
be constrained by photometry if the BH mass and accre- 
tion rate can be independently estimated (e.g. Czerny et al. 
2011). This method so-far has large uncertainty mainly due 
to the large uncertainty of BH mass. 

The comparison between a quasar's UV emission and 
its polarized visible emission can be complicated by the 
possible Balmer edge absorption. It has been reported that 
the polarized flux of the two quasars manifest a disconti- 
nuity in the slope at the wavelengths shorter than ~400 
nm (Kishimoto et al. 2004). The feature is mostly inter- 
preted as the buried Balmer edge of the intrinsic spectra of 
the quasar UV/optical continuum, or the "Big Blue Bump" 
(BBB) emission. The edge absorption feature is due to the 
bound-free opacities in the disc atmosphere, and indeed in- 
dicates the thermal and optically thick nature of the contin- 
uum. However, the origin of the Balmer edge absorption is 
still unknown. An order-of-magnitude estimate of the broad- 
ening of the Balmer edge indicates that the broadening is 
consistent with the orbital motion of the corresponding disc 
annuli responsible for emission around 400 nm (Kishimoto 
et al. 2004). 

In this work, we use the extended standard thin ac- 
cretion disc model to fit the polarized spectra in the vis- 
ible and near-infrared bands and the total spectra in the 
UV band of 3C 95 and 4C 09.72, after removing prominent 
emission and absorption features; the extension includes the 
Balmer edge absorption of an optically thin layer of neutral 
gas co-rotating with the accretion disc. We use the Monte- 
Carlo-Markov-Chain (MCMC) method to fully explore the 
multi-dimensional parameter space of BH mass, spin, accre- 
tion rate, disc inclination and polarization. The wide spec- 
tral coverage of the total UV emission and the polarized 
spectra in visible and NIR bands provides a complete de- 
scription of the accretion disc emission, which allows us to 
determine of physical properties of the BH and the accretion 
disc. In particular, the true polarization fraction of the disc 
emission can be well constrained from the spectral fitting, 
which leads to a continuum spectrum from UV to NIR. The 
spectral fitting also constrains the broadening of the Balmer 
edge, which is constroled by the disc inclination. 
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2 STANDARD ACCRETION DISC MODEL 
WITH BALMER EDGE 

We intend to describe the polarized spectra in visible and 
NIR wavelengths as well as the total spectra in UV wave- 
lengths by the standard thin disc model (Shakura & Sunyaev 
1973; Lynden-Bell & Pringle 1974). As an extension to the 
standard model, we include the Balmer edge absorption of 
the neutral gas that is assumed to be co-rotating with the 
accretion disc. We do not intend to model in detail the atmo- 
sphere of quasars' accretion discs (e.g., Loar & Netzer 1989; 
Loar et al. 1990; Hubeny et al. 2000) due to the extensive 
computation burden; our main goal is to determine the BH 
parameters from quasars' continuum. 

According to the standard thin disc model, its temper- 
ature of the accretion disc follows (Krolik 1998; Hubeny et 
al. 2000) 
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in which T is the temperature, G is the gravitational con- 
stant, A/bh is the BH mass, M is the accretion rate, a is 
the Stefan-Boltzmann constant, r is the radius, x is the ra- 
dius in the unit of gravitational radius r g = ga/ bh , r\ is the 
radiative efficiency as function of BH spin a/M, L/Ledd is 
the accretion rate in the unit of Eddington accretion rate, 
and J?r is the correction term of GR effects analytically ex- 
pressed as a function of a/M (see Krolik 1998). Analytical 
expressions of rj and _Rr are given in Appendix (fXJ). 

We assume that the Balmer edge absorption comes from 
a relatively narrow disc annuli for effective temperature from 
8,000 to 20,000 K, based on the consideration of the n = 2 
population of hydrogen. The ionization of hydrogen depends 
on the number density in the photosphere, which in turn de- 
pends on the torque of the disc, an observational unknown. 
With a sensible range of the torque value, a realistic disc 
vertical structure model (Laor & Netzer 1989), and the lo- 
cal thermal equilibrium of a pure hydrogen plasma, we find 
that the n — 2 population concentrates in a relatively nar- 
row range whose temperature is from 8,000 to 20,000 K. A 
large electron density from ionization of metals, and a range 
of other parameters might affect the temperature range of 
n = 2 population; we have performed a sensitivity study 
by letting the Balmer edge absorption to occur in annuli of 
T = 8, 000 ~ 15, 000 - 25, 000 K. The frequency-dependent 
form of the Balmer edge feature is, 
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where to is the characteristic optical depth due to the com- 
mon photoelectric absorption, v c is the threshold frequency 
corresponding to ~ 364.6 nm, and /i = cos(i) in which i is 
the inclination angle of the accretion disc. We assume to to 
be a constant in the Balmer absorption annuli and zero else- 
where for simplicity. Given L, Mbh, and a/M, the range of 
radii of Balmer absorption annuli xo can be found by letting 
T = 8, 000 ~ 20, 000 K in Eq. (Q. 

The inclination of the accretion disc has three effects: 
local limb darkening, projection of the disc on the sky plane, 



and the relativistic doppler effect due to the orbital motion. 
Omitting the general relativistic effect of photon propaga- 
tion, the asymmetric broadening factor due to the orbital 
motion is derived as 
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where v is the observed frequency, vo is the emission 
frequency in the rest frame of the rotating gas, 70 = 
1 / \J 1 — jSq , P\ is the relativistic factor corrected for the in- 
clination as /3i = /?osin(i), /3o = vo/c, and vo is the local 
Keplerian velocity of the disc where the Balmer absorption 
takes place. Here we see that the broadening of the Balmer 
edge feature sensitively depends on the disc's inclination. In 
Newtonian regime, /9 = x -1 ^ 2 . In the Kerr spacetime, is an 
analytical function of x and a/M (see Krolik 1998 for detail 
formulation). We verify that when x is large, the functional 
form of f3 in the Kerr spacetime approaches to that of the 
Newtonian regime. 

The emerging spectrum from an optically thick disc 
with an optically thin layer of absorbing neutral gas is then 
described by 
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where the local limb darkening and projection of the disc 
plane are considered as in Krolik (1998). The inner edge of 
the accretion disc (xi n ) is assumed to be located at the inner 
most circular orbit of the spinning BH, as a function of a/M. 
The analytical expression of Xi n is given in Appendix ((A)) . 
In our calculation we choose a; ou t = 1000, large enough to 
compute the NIR spectra up to 2 /jm. f(x, v), the emission 
flux from an annuli between x and x + dx, can be computed 
from the following integration 

f(x,v) = / B(T(x), v ) exp(-r(v ))S(v, v )dv , 

(5) 

where B(T, v) is the Planck function. The observed polar- 
ized spectrum in the optical and NIR wavelengths {vF Vt om) 
and the UV continuum (yFu,w) are then 

vF„,om = P-jjT . "^uv = , (6) 

where P is the polarization fraction, and D is the luminosity 
distance. The luminosity distance to a source is calculated 
from the redshift z with the standard flat cosmology model 
of (Om = 0.3, = 0.7) and Hubble constant Ho ~ 72 km 
s" 1 Mpc" 1 . 

In summary, our model computes the continuum spec- 
trum from UV to NIR of the thin accretion disc with the 
effect of co-rotating optically thin Blamer edge absorption. 
The model involves 6 independent parameters: BH mass 
(Mbh), BH spin (a/M), accretion rate (L), polarization (P), 
inclination (i), and optical depth of the Balmer absorption 
(to). In this work we use the kinetics of BLR and the veloc- 
ity dispersion to determine the BH mass, which decreases 
the number of free parameters to 5. We further define the 
"effective" accretion rate as L' = L x + |)ju according 
to Equation Q, in order to remove the correlation between 
the accretion rate and the inclination. Finally, a/M, L' , P, i 



4 R. Hu & S. N. Zhang 



and ro are our five independent fitting parameters. We em- 
phasize that the formulation here is a simplified treatment 
of quasars' accretion disc emission continuum. This formu- 
lation is computationally efficient, which allows extensive 
exploration of the parameter space with the Monte-Carlo- 
Markov-Chain method. 



3 DATA SELECTION 

We obtain the total-light and polarized spectra of Q0144- 
3938, 3C 95, CTS A09.36, 4C 09.72, PKS 2310-322 and Ton 
202 in the visible and NIR wavelengths from Kishimoto et 
al. (2008). The optical spectropolarimetry was taken with 
the instrument FORS1 mounted on the Very Large Tele- 
scope (VLT) UT2 and the Low Resolution Imaging Spectro- 
graph (LRIS) on the Keck-I telescope (Kishimoto et al. 2004; 
Kishimoto et al. 2008). 3C 95 and 4C 09.72 were observed 
with the VLT in September 2002. Q0144-3938, CTS A09.36, 
PKS 2310-322 were observed with the VLT in September 
2005. Ton 202 was observed with the Low Resolution Imag- 
ing Spectrograph (LRIS) on the Keck-I telescope on May 4, 
2003. 

The NIR broad-band imaging polarimetry of the 
six quasars was obtained with the instruments UFTI 
(with the polarimetry module IRPOL2) and ISAAC which 
are mounted on the United Kingdom Infrared Telescope 
(UKIRT) and VLT UT1, respectively. Ton 202 was observed 
with the UKIRT on January 16, 2001 (Kishimoto et al. 
2005). 3C 95 and 4C 09.72 were observed with the UKIRT 
in the fall of 2006. Q0144-3938, CTS A09.36, PKS 2310-322 
were observed with the VLT on August 21, 2007. 

Among the sample of the six quasars, 3C 95 and 
4C 09.72 have reduced spectra in the near UV wavelengths 
available for scientific interpretations. The total-light spec- 
tra of 3C 95 and 4C 09.72 in the UV wavelength were taken 
by the HST Faint Object Spectrograph (FOS) and avail- 
able from the archive of the HST High-Level Science Prod- 
ucts (Evans & Koratkar 2004). The spectrum of 3C 95 was 
taken on December 21, 1991 with a total exposure time of 
4454 s. The spectrum of 4C 09.72 was taken on October 12, 

1992 with a total exposure time of 4751 s and on July 7, 

1993 with an exposure time of 480 s, respectively. The UV 
and optical spectra are recalibrated uniformly with all other 
HST pres-COSTAR (Corrective Optics Space Telescope Ax- 
ial Replacement) FOS spectrophotometric data for AGNs 
and quasars, using the up-to-date algorithms and calibra- 
tion data (Evans & Koratkar 2004). Where possible, multi- 
ple observations are combined to produce a single spectrum 
for each object with the highest possible signal-to-noise ra- 
tio and covering the widest wavelength range. We retrieve 
the recalibrated spectra of 3C 95 and 4C 09.72 in 1140-3300 
A in this work. 

The retrieved UV spectra have many prominent emis- 
sion and absorption lines, and have the spectral resolution of 
~1 A. The UV spectra have been processed in the following 
steps before taken into the spectral fittings. First, we remove 
the absorption lines in the UV spectra. For 3C 95, we use 
the complete sample of absorption lines provided by Bah- 
call et al. (1993), which contains the central wavelength in 
the observation frame and FWHM of each absorption line. 
We assume a Gaussian profile for each absorption line with 



Table 1. Absorption lines removed from the UV spectra of 3C 95. 
The wavelength is in the observation frame. 



A 


FWHM 


\F\ x 10 11 


Line 


(ran) 


(nm) 


erg s _1 cm~ 2 ) 




1649.02 


1.51 


0.30 


Lyman a 


1742.76 


1.86 


0.18 


Lyman a 


1749.42 


1.65 


0.32 


Lyman a 


1890.00 


2.08 


0.35 


Si IV 


1909.82 


1.88 


0.56 


Lyman a 


1959.37 


1.59 


0.58 


Lyman a 


2025.72 


1.94 


0.25 


Zn II (Mg II) 


2099.70 


1.84 


0.53 


C IV 


2103.38 


1.80 


0.43 


C IV 


2314.49 


4.64 


0.14 




2382.53 


2.04 


0.29 


Fe II 


2586.82 


2.04 


0.18 


Fe II 


2600.20 


2.13 


0.31 


Fe II 


2755.34 


2.04 


0.17 




2796.25 


2.23 


0.35 


Mg II 


2803.27 


2.55 


0.25 


Mg II 


2852.95 


2.04 


0.19 


Mgl 



three parameters: central wavelength, FWHM, and strength. 
We take the central wavelength and the FWHM from Bah- 
call et al. (1993), and fit the UV spectrum to determine the 
strength. We list in Table [1] central wavelengths, FWHM, 
and strengths of all absorption lines. For 4C 09.72, because 
no absorption line list is available from the literature, we 
do not perform the removal of absorption lines; because the 
absorption lines are all very narrow with FWHM less than 
3 A, they are unlikely to affect the fitting to the continuum. 

Second, we remove the emission lines in the UV spec- 
tra. The extraction of emission lines, especially the promi- 
nent broad lines, are complicated. Marziani et al. (1996) dis- 
cussed the profiles of emission lines of 3C 95 and 4C 09.72 
in details. In this work, however, we only focus on the con- 
tinuum. Hence, we try to remove the emission lines in a 
simple way. We use the line identifications and central wave- 
lengths given by the HST High-Level Science Products and 
remove each emission line by several trial Gaussian profiles. 
Note that there is no spectral fitting performed here. We 
test each of the emission lines with several combinations of 
FWHM and line strength, assuming the Gaussian profiles, 
and determine the best line removal by manual inspection. 
The Gaussian profiles retracted from the UV spectra are 
tabulated in Table [2] We do not attempt to remove the iron 
emission lines. For both 3C 95 and 4C 09.72, the Fe II and Fe 
III pseudocontinua are much weaker than the C IV emission 
(Marziani et al. 1996). Also, we only use the UV continua 
for the rest-frame wavelengths shorter than ~ 200 nm; in 
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Table 2. Emission lines removed from the UV spectra of 3C 95 
and 4C 09.72. The wavelength is in the observation frame. 



Object 


A 


FWHM 


\F X x 10 11 


Line 




(nm) 


(nm) 


(erg s~ 1 cm -2 ) 




3C 95 


1670.9 


25 


0.30 


Ly/3, O VI 


3C 95 


1964.9 


70 


0.80 


hya 


3C 95 


1964.9 


30 


0.50 


Lya 


3C 95 


2004.5 


60 


0.30 


N V 


3C 95 


2257.4 


50 


0.10 


Si IV 


3C 95 


2504.4 


90 


0.50 


C IV 


3C 95 


2504.4 


35 


0.30 


C IV 


3C 95 


3084.9 


80 


0.20 


C III] 


4C 09.72 


2001.9 


40 


0.15 


Si IV 


4C 09.72 


2220.9 


80 


0.50 


C IV 


4C 09.72 


2220.9 


20 


0.75 


C IV 


4C 09.72 


2735.7 


40 


0.20 


C III] 



this wavelength range the effect of iron lines to our spectral 
fitting is minor (e.g., Vestergaard & Wilkes 2001). Finally, 
for 3C 95 and 4C 09.72, we dismiss the spectrum with wave- 
length shorter than 1300 A in the rest frame because of 
the prominent Galactic Lyman-a emission and large statis- 
tical errors (e.g. Evans & Koratkar 2004), the non-thermal 
component due to Compton scattering in the accretion disc 
atmosphere (e.g., Hubeny et al. 2001), and possible inter- 
galactic contamination (Telfer et al. 2002). The UV contin- 
uum spectra after the absorption and emission line removal 
are shown in Figure [T] 

Finally, we transform the spectra into the rest frame by 
multiplying the flux F\ by (1 + z) and present the flux in 
terms of \F\ (see blue curves in Figure [TJ. For the fitting 
purpose we also bin the UV spectra with a linear wavelength 
spacing of 1 nm, as shown in Figure [T] The uncertainties of 
the flux are combined by the standard error propagation 
method. 



4 DATA ANALYSIS 

4.1 Black Hole Mass Determination 

The most direct measurements of black masses come from 
the kinematics of the broad line regions (BLRs) of AGNs, 
i.e., the emission-line width SV and the size of BLR given by 
reverberation mapping (Blandford & Mckee 1982; Peterson 
1993), assuming that the line emitting gas in the BLR is 
varialized. The virial mass can then be expressed as 



where / is a factor depending on the structure, kinemat- 
ics and orientation of the BLR (Peterson et al. 2004). On 
average, the factor / can be taken as < / >= 5.5 if the 
line width is taken as the line dispersion (Time- For a distant 
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Figure 1. Pre-fitting processes of the UV spectra of 3C 95 and 
4C 09.72. Blue curves show the original spectra retrieved from the 
archive of the HST High-Level Science Products with wavelengths 
in the rest frame. Red curves show the spectra after the emission 
and absorption lines removed. Black curves show the spectra with 
a linear wavelength spacing of 1 nm. 



AGN the reverberation measurement is not possible; then 
the BLR size comes from the radius-luminosity (R-L) rela- 
tionship (the relationship between the radius of BLR and the 
optical luminosity) calibrated at low redshift (Kaspi et al. 
2005) . The removal of the host-galaxy starlight contribution 
improves the relationship significantly, and the relationship 
manifests as r oc L 1 ^ 2 (Bentz et al. 2006). 

The BH mass determined this way has been found to be 
tightly correlated with the velocity dispersion of the bulge 
or spheroid <r» (the so-called Mbh — <j* relation; Ferrarese & 
Merritt 2000; Tremaine et al. 2002). In this work, the Mbh — 
a* relation we used is determined from the reverberation- 
mapped virial BH mass, which is provided by Onken et al. 
(2004). In general, the FWHM of the [O III] A5007 emission 
line can be used as a proxy for a, (Nelson & Whittle 1996; 
Nelson 2000). 

Although we only perform the spectral fitting to the 
continuums of 3C 95 and 4C 09.72 in the following, we 
determine the BH masses for all six quasars in the sam- 
ple from the BLR kinematics and compare them with the 
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10.0 




Figure 2. The virial black hole masses Mjh and the velocity 
dispersion <r* on A/bh — c* relations. The solid line shows the 
A/bh — c» relation by Tremaine et al. (2002). The uncertainty 
of Mbh is within a factor of 2 (Peterson et al. 2004), and the 
uncertainty of log(<r„) is 0.07 dex (Tremaine et al. 2002). 



Mbh — c* relation. In this way we can verify the reliability 
of the BH mass determination. We obtain the broad H/3 line 
width by interpolating a linear continuum between contin- 
uum windows on either side of the line, and calculate the 
line's mean squared dispersion defined as the second mo- 
ment of the line profile (Peterson et al. 2004). In practice, 
the narrow-line residuals influence significantly the deter- 
mination of the broad-line width. Therefore we remove the 
narrow-line component carefully using the adjacent [O III] 
A5007 line as a template. The obtained line widths and BH 
masses are tabulated in column 1 of Tabled The BH masses 
determined by such method have typical errors by a factor 
of 2, i.e., 0.3 dex (Peterson et al. 2004). 

Using the line dispersion of [O III] A5007 as the probe of 
the velocity dispersion of an AGN's host galaxy's bulge, we 
can plot the obtained BH masses with respect to the estab- 
lished Mbh — cr* relation (Tremaine et al. 2002), as in Figure 
[5] Using the Nukers' estimate, we find the reduced x 2 is 1-2, 
if errors for A/bh-rm and er„ are taken as 0.3 dex (Peterson 
et al. 2004) and 0.07 dex (Tremaine et al. 2002), respectively. 
Mbh-rm measures the virial mass the BH, and <t» indicates 
the dispersion mass of BH. The consistence of these inde- 
pendent measurements as shown in Figure [2] allows us to 
conclude that BH mass estimates with both methods are 
reliable for these sources. To further improve the accuracy 
of these BH's mass estimates, we combine the BH mass es- 
timates from both methods; the errors for the combined BH 
masses are thus 0.21 dex. In Table 1, we list the observation 
data used and the BH masses determined from the reverber- 
ation mapping (denoted as Mbh-rm), the velocity disper- 
sion (denoted as A/bh-<t), and the combined kinematic mea- 
surements log Mbh-km = 0.5 x (log Mbh-rm +log A/ B h-<t)- 



4.2 Spectral fitting 

We fit the two observed quasar continuum spectra (3C 95, 
4C 09.72) to our extended standard think disc model pre- 
sented above, with the polarized spectra in the visible and 
NIR wavelengths, and the total spectra in the UV wave- 
lengths. The polarized spectra are bridged with the total 
light spectra by the polarization fraction of the observed ac- 
cretion disc emission. The 5 fitting parameters are BH spin 
o/M, accretion rate L' , polarization fraction P inclination 
i, and Balmer edge optical depth to. In the fitting, we use 
the BH mass estimate Mbh-km in table [3] 

We use the Monte Carlo Markov Chain (MCMC) 
method for the spectral fitting to explore the whole param- 
eter space. Due to the multi-parameter nature of the prob- 
lem, the MCMC method not only provides possibility of fast 
computation, but also gives the posterior probability distri- 
bution of all parameters. We use the adaptive Metropolis- 
Hastings algorithms implemented by Haario et al. (2006). 

The convergence of the Markov chain is examined as 
follows. We start from two randomly selected parameter sets 
and generate two Markov chains. The second halves of both 
chains are mixed together and the so-called "potential scale 
reduction" (7?) is computed for each estimated parameter 
(Gelman & Rubin 1992); R declines to unity as the length 
of Markov chain goes to infinity. In this work we require R < 
1.005 for all estimated parameters as the criterium for the 
convergence of Markov chain. We find that this requirement 
corresponds to a Markov chain length of 1 millions, so for 
each quasar at least 2 millions continuum spectra have been 
generated to be compared with the data. 

The fitting results are tabulated in Tableland the best 
fitted model spectra are shown in Figure [3] as solid lines. As 
comparison, we also plot the corresponding model spectra 
without the Balmer-edge absorption as dashed lines. The 
probability density of each parameter from the MCMC sim- 
ulation is shown in Figuref?] The optimal value (Xo) of each 
parameter is taken as the best-fitted value of the parameters 
(i.e., the parameter set that provides the smallest \ 2 m the 
Markov chain). Due to the asymmetry of these probability 
distributions, the lower (Al) and upper (Au) boundaries at 
95% confidence level of each parameter is determined in the 
following way: 

r x o r x o 

/ p(X)dX/ / p(X)dX = 0.95, (8) 
'x L Jo 

and 

rX v roo 

/ p(X)dX/ / p(X)dX = 0.95, (9) 

J X Q J Xq 

where p(X) is the probability density of the parameter A. 
The correlations among the five fitting parameters of 3C 95 
and 4C 09.72 are shown in Figure [5] and Figure [(J respec- 
tively. We have provided both the best fitted values and the 
most probably values for the parameters in Table 21 The 
discrepancies between the best-fitted values and the most- 
probable values can be quite large, for those parameters that 
are poorly constrained by the data (e.g., a/M, i). Neverthe- 
less, for all parameters the best-fitted values and the most- 
probable values lie within their 1-a ranges. 
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Table 3. Quasar black hole mass estimates 



Quasar name 


H/3 width 
(km s — 1 ) 


BLR size 1 
(Light days) 


Mbh-rm 2 
(10 8 M Q ) 


<7, 3 

(km s" 1 ) 


Mbh- ct 4 
(1O 8 M ) 


A/bh-km 5 
(10 8 M Q ) 


z 


Q0144-3938 


2570 


55.8 


3.94 


362 


14.7 


7.6 


0.244 


3C 95 


2846 


171.5 


14.9 


394 


20.6 


17.5 


0.616 


CTS A09.36 


3124 


54.8 


5.72 


218 


1.9 


3.3 


0.310 


4C 09.72 


2855 


179.8 


15.7 


359 


14.2 


14.9 


0.433 


PKS 2810-322 


3048 


70.4 


7.00 


197 


1.3 


3.0 


0.337 


Ton 202 


3306 


77.4 


9.05 


380 


17.8 


12.7 


0.366 



1 The r-L relationship of Kaspi et al. (2005) is applied. 

2 Estimated from H/3 width and BLR size, with a 1-cr error of 0.3 dec. The factor / is taken as / = 5.5. 

3 The velocity dispersion is taken as the line dispersion of [O III] A5007. 

4 Estimated from cr t , with a 1-<t error of 0.3 dec. 

5 Geometrical average between A/bh-rm an d A/bh- CT i i.e., log A/bh— KM = 0.5 X (logA/sH— rm + l°g A/bh— <t)> with a 1-cr error of 
0.21 dec. 

Table 4. Physical parameters of 3C 95 and 4C 09.72 estimated by fitting their continuum spectra. Five parameters are estimated from 
the spectral fitting: black hole spin a/M, accretion rate L'/^EDD = i/-^EDD x §(a* + §)^i polarization P, disc inclination i and Balmer 
edge optical depth to. The spectral fitting is performed with MCMC simulations that generate chains to the convergence. The black hole 
masses are estimated by kinetic measurements (A/bh— km in Table [3}. The first line of each object lists the best fitted results, defined 
as the parameter set that provides the smallest x 2 . The second line of each object lists the most probable parameter values, defined as 
the peak value in the posterior probability distribution. The third line of each object lists the lower and upper 95% boundary values of 
parameter estimates on both sides of the best-fitted values for T = 8,000 ~ 20,000 K, using Eqs. (8) and (9), respectively. The fourth 
line lists the parameter boundary values estimated from a separate set of MCMC simulations, in which we have assumed the Balmer edge 
absorption occurs in the disk's annuli of T = 8, 000 ~ 15, 000 K. Similary, the fifth line lists the parameter boundary values estimated 
with T = 8,000 ~ 25,000 K. We see that except for tq, the spectral fitting does not depend on the specific choice of the Balmer edge 
temperature range. 





a/M 


i'/^EDD 


P (%) 


i (°) 


TO 


X 2 /dof 




3C 95 (Best-fitted) 


0.030 


0.070 


1.48 


70.4 


0.32 


4237.8/2637 


0.28 


Peak 


0.65 


0.067 


1.37 


51.6 


0.45 




0.11 


T = (8 — 20) x 10 3 K 


0.0035 - 0.84 


0.061 - 0.080 


1.08 - 1.83 


6.6 - 79.6 


0.035 - 4.80 




0.049 - 0.70 


T = (8 - 15) X 10 3 K 


0.0054 - 0.82 


0.062 - 0.084 


1.10 - 1.86 


8.1 - 82.7 


0.055 - 9.51 




0.050 - 1.11 


T = (8 - 25) x 10 3 K 


0.0045 - 0.84 


0.061 - 0.081 


1.10 - 1.82 


6.8 - 80.0 


0.021 - 1.60 




0.050 - 0.75 


4C 09.72 (Best-fitted) 


0.0049 


0.021 


2.46 


67.1 


0.34 


7678.7/3295 


0.067 


Peak 


0.48 


0.020 


2.09 


30.4 


0.75 




0.020 


T = (8 — 20) x 10 3 K 


0.0024 - 0.76 


0.018 - 0.027 


1.49 - 3.25 


3.5 - 78.3 


0.035 - 4.04 




0.014 - 0.20 


T = (8 — 15) x 10 3 K 


0.0034 - 0.74 


0.018 - 0.024 


1.60 - 3.12 


3.9 - 75.3 


0.135 - 9.56 




0.015 - 0.14 


T = (8 - 25) x 10 3 K 


0.0034 - 0.77 


0.018 - 0.027 


1.47 - 3.09 


3.3 - 74.8 


0.036 - 1.96 




0.015 - 0.15 



5 RESULTS 

5.1 Polarization of Accretion Disc Emission 

The continuum spectra of 3C 95 and 4C 09.72 from UV to 
NIR wavelengths can be sufficiently described by the stan- 
dard accretion disc model with a prominent Balmer-edge 
absorption, as shown in Figure [3] The peak emissions of 
both 3C 95 and 4C 09.72 are at the wavelengths of 100~200 
nm, which have been observed by the HST. The observed 
UV continuum spectra of the two quasars are indeed the 
optically thick thermal radiation from the inner part of the 
accretion disc. The polarized spectra from visible to NIR 
wavelengths are wavelength-independently scaled spectra of 
the disc's thermal emission; therefore the polarized spectra 
from visible to NIR wavelengths can be bridged with the to- 
tal UV continua by adjusting the polarization fraction. The 



resulting continuum spectra from UV to NIR represent the 
thermal emission from a wide radius range of the quasars' 
accretion discs. 

Optical polarization of accretion disc emission can be 
determined by bridging the polarized visible-IR emission 
to the total emission at UV wavelengths. Both total emis- 
sion and polarized emission, although at different wave- 
lengths, follow the same multi-color thermal emission spec- 
trum; hence the polarization fraction can be well constrained 
by spectral fitting of the continuum thermal emission of ac- 
cretion discs (see Figure [4} . The uncertainty of the opti- 
cal polarization is gaussian and independent from other BH 
or accretion disc properties (see Figure [5] and Figure [6}. 
We emphasize that the polarization derived in this paper 
is the true polarization of the accretion disc emission, and 
we expect the true polarization to be larger than the ob- 
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X [nm] 



400 600 800 1000 
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Figure 3. Best fitted model spectra of 3C 95 and 4C 09.72. The upper panels show the observed spectra and the fitting results. The 
observed polarized spectra have been uniformly multiplied by a factor of l/P to be compared with the UV total light continuum. The 
solid lines are the best fitted model spectra with the parameters tabulated in Table [4] The dashed lines are calculated with the same 
parameters but without the Balmer-edge absorption (to = 0). The lower panels show the fitting residuals. Note that the vertical axis of 
the lower panels is in the linear scale. We see that the Balmer edge absorption is clearly revealed in the quasars' polarized continuum 
spectra. 



served polarization due to host galaxy contamination. The 
observed polarization fraction of both 3C 95 and 4C 09.72 
is ~ 1.0 — 1.5% at wavelengths of 400 ~ 500 nm and gradu- 
ally decreases to ~ 0.5 — 1.0% at 250 nm (Schmidt & Smith 
2000; Kishimoto et al. 2008). These observations are con- 
sistent with our spectral fitting results, although we do not 
use total light in visible wavelengths directly. In comparison 
with the true polarization obtained from spectral fitting, we 
find that for 3C 95, there is little contamination from its 
host galaxy at 400 ~ 500 nm, but the contamination in- 
creases significantly for shorter wavelengths; for 4C 09.72, 
host galaxy contamination is significant for all visible and 
NIR wavelengths. 

Spectral fitting to the accretion disc continuum from 
UV to NIR may constrain the BH and accretion disc prop- 
erties. The disc's accretion rate is tightly constrained by 
the total luminosity of the accretion disc, a physical quan- 
tity that is derived from the polarized luminosity and the 
true polarization fraction (see Figure [4]). The accretion rate 
is only weakly correlated with the BH spin, as predicted 



by the standard thin disc model (see Figure [5] and Figure 
Also, the estimated accretion rates are significantly sub- 
Eddington, well within the range that the standard thin disc 
model can be applied reliably. We have also performed con- 
tinuum spectral fitting for 3C 95 and 4C 09.72 with the BH 
mass assumed as the fitting parameter (not shown in the 
paper). We obtained quantitatively very similar results and 
the BH mass estimated from the spectral fitting is consis- 
tent with the kinetics estimates in Table [3] This suggests 
that their disc continuum spectra do not contain additional 
information on their BH masses, i.e., their BH mass uncer- 
tainties are dominated by the uncertainties in their kine- 
matic mass estimates. 



5.2 Revealed Balmer Edge Absorption 

Prominent Balmer edge absorption features are found in 
their continuum spectra, as shown in Figure [3] The Balmer 
edge absorption is due to the bound-free opacities in the 
accretion disc atmosphere, which can only be seen in the 
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Figure 4. Posterior probability densities of physical parameters of 3C 95 and 4C 09.72. For each object, the probability densities 
are computed from the MCMC chain of the last 2 millions iterations to minimize the effect of initial conditions. The accretion rate 
(L/Ledd) is computed from the effective accretion rate (Z/) and the inclination (i) of the fitting results. We see that UV continuum 
and polarized optical continuum can well constrain the accretion rate and the polarization, and the balmer edges revealed by polarized 
optical continuum provides information of the inclination. 
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disc's emission continuum with host galaxy contamination 
removed. We here confirm that the Balmer edge absorption 
for quasars is observable via polarized spectra in the visible 
wavelengths (i.e., Kishimoto et al. 2004; 2005; 2008). The 
predicted shape of the Balmer edge absorption combines the 
photoelectric absorption profile and the broadening due to 
the orbital motion. As shown in Figure [3j the Balmer edge 
absorption fits the break of polarized spectra at ~ 400 nm 
and connects to the UV continuum smoothly. 

The observed Balmer edge absorption is broadened due 
to the orbital motion of the accretion disc. Since the Balmer 
edge opacity is proportional to the concentration of n = 2 
hydrogen, the Balmer edge absorption feature in the disc 
continuum is generated in a narrow annuli of the accre- 
tion disc in which the temperature is 8, 000 ~ 20, 000 K. 
The inclination of the accretion disc affects significantly the 
broadening of the Balmer edge feature, as the broadening 
mostly depends on the line-of-sight projection of the orbital 
motion velocity of the specific Balmer-edge-generating an- 
nulus. From the spectral fitting, we find that for one thing, 
orbital motion is sufficient to account for the broadening of 
the revealed Balmer edges for 3C 95 and 4C 09.72 (Figure 
[5J; for another, the widths of Balmer edges may constrain 
the discs' inclination (Figure |4)|. Noticeably, the edge-on ge- 
ometry is confidently ruled out by fitting emission spectra 
for 3C 95 and 4C 09.72 (FiguregJ, consistent with them be- 
ing type-1 quasars. However, the spectral fitting cannot rule 
out the face-on geometry (Figure [4}, because of the intrinsic 
width of the Balmer edge from the v/u~ 2 ' 67 dependency of 
the photoelectric absorption cross sections. One might as- 
sert that the most-likely inclination of 3C 95 is larger than 
that of 4C 09.72; any further constraint of disc inclination 
by Balmer edge broadening is impeded by a relatively low 
signal-to-noise ratio of the Balmer edge feature, and the un- 
certainties of Balmer edge optical depths. 

Spectral fitting to the accretion disc continuum from 
UV to NIR can determine the optical depth of the Balmer 
edge absorption. As shown in Tabled the Balmer edge opti- 
cal depths of 3C 95 and 4C 09.72 are large enough to create 
significant spectral features. The upper limit of the optical 
depth is constrained by the spectral fitting, indicated by 
the posterior probability distribution shown in Figure [4] We 
caution that for the Balmer edge optical depth significantly 
larger than unity the optically thin assumption may break 
down and the full radiative-transfer model of the accretion 
disc atmosphere is needed. After all, the best-fitted values 
and the most-probable values of the optical depth for both 
quasars are found to be consistent with an optically thin 
absorption. With the bound-free cross section of hydrogen 
at the Balmer edge (ffbf = 1.3 x 10~ 17 cm 2 ), we estimate 
the column density of excited (n = 2) hydrogen atoms to be 
3x 10 15 - 4x 10 17 cm" 2 for 3C 95 and 3x 10 15 - 3x 10 17 cm" 2 
for 4C 09.72, respectively. For local thermal equilibrium the 
implied neutral hydrogen column density on the accretion 
disc by the Balmer edge optical depth is estimated to be 
10 19 ~ 10 21 cm . The neutral hydrogen column density is 
comparable to the neutral hydrogen column density of the 
Milky Way, which is consistent with the interpretation that 
the Balmer edge features originate on the quasars' accretion 
discs. 



6 CONCLUSION AND DISCUSSION 

A quasar's polarized continuum spectrum at visible to NIR 
wavelengths, together with the total continuum spectrum 
at UV wavelengths, provides an unique opportunity to de- 
termine the properties of accretion discs around supermas- 
sive BHs. In this work we fitted the spectra of two type-1 
quasars (3C 95 and 4C 09.72) to the standard thin accretion 
disc model with the Balmer-edge absorption, with the BH 
masses determined from stellar dynamics and reverberation 
measurements. We used the MCMC method to sample a pa- 
rameter space of BH spin, accretion rate, polarization, disc 
inclination and Balmer edge optical depth. 

We conclude that the quasar polarized spectra at visible 
and NIR wavelengths and total spectra in UV wavelengths 
can reveal the accretion disc emission, which is indeed the 
thermal emission of the optically thick and geometrically 
thin accretion disc. At the UV wavelengths, the contamina- 
tion from the host galaxy is negligible so the total spectra 
reflects the emission spectra of the accretion disc. At the 
visible and NIR wavelengths, the contamination buries the 
intrinsic accretion disc emission as well as the Balmer-edge 
feature. However the contamination does not exist in the 
polarized spectra. Therefore, the polarized spectra at vis- 
ible and NIR wavelengths and total spectra at UV wave- 
lengths allow us to study the entire accretion disc. The true 
optical polarization of the accretion disc emission can be 
well constrained by bridging the total emission at UV and 
the polarized emission at visible. We find that 3C 95 and 
4C 09.72 have significantly sub-Eddington accretion rates, 
and their optical polarization is 1.1 — 1.8% and 1.5 — 3.3% (at 
95% confidence level), respectively. The true optical polar- 
ization of accretion disc derived from our continuum fitting 
is larger than the observed optical polarization, which in- 
dicates host galaxies' contamination in the quasars' total 
optical continua. 

The absorption feature found at the Balmer edge is 
prominent and significantly broadened. We found that the 
Balmer absorption feature for both quasars investigated can 
be reproduced by assuming optically thin absorbing n — 2 
hydrogen population co-rotating with the accretion disc at 
radii with effective temperature between 8,000 and 20,000 
K. The Balmer edge feature is found to extend to UV wave- 
lengths, whose width is the result of both photoelectric ab- 
sorption cross section and the line-of-sight projection of the 
orbital motion. The choice of the temperature range does 
not significantly affect the spectral fitting. For example, for 
the upper boundary of temperature varying from 15,000 to 
25,000 K, we did not find any significant corresponding de- 
pendency of the black hole parameters and the true optical 
polarization (see Table|3|. However, we found the best-fitted 
Balmer edge optical depth (to) to be strongly dependent on 
the width of Balmer absorption annulus. The wider the an- 
nulus is, the deeper the Balmer edge appears, and the re- 
quired to to fit the data is smaller (see Table [4}. Therefore 
the estimation of neutral hydrogen on the surface of the ac- 
cretion disc is also sensitive to the temperature range. The 
sensitivity study has suggested that a reasonable uncertainty 
in the exact location where the Balmer edge absorption oc- 
curs has a major effect in the depth of the Balmer edge fea- 
ture, but only has a minor effect on the broadening of the 
Balmer edge feature. Indeed, the broadening of the Balmer 
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Figure 5. Correlation of the five fitting parameters of 3C 95. Each panel shows the joint posterior probability distribution between a 
pair of parameters with each dot corresponding to a parameter set in the MCMC chain of the second half of the 2 million iterations. 



edge feature depends mostly on the disc's inclination. The 
broadening of the Balmer edge features tentatively constrain 
the disc's inclination of the two quasars. 

Our fitting results show that the quasars' BH spin pa- 
rameters cannot be constrained well from the spectral fitting 
up to ~ 100 nm. As shown in Figure [4] and Table O the pos- 
terior probability distribution of the BH spin is very wide. 
The best-fitted parameters for both 3C 95 and 4C 09.72 in- 
dicate non-rotating BHs; however, the full exploration to the 
parameter space indicates that a significant range of a/M 
between to 1 is acceptable. This is due to the combination 
of the uncertainties in inclination, accretion rate and BH 
mass. 

Figure [S] and [5] show a weak correlation between to and 
i. As expected, a larger to would produce a deeper absorp- 
tion edge, which allows the widths of the Balmer edge fea- 
ture to be better determined, and the disc's inclination to be 
better constrained. Indeed, the uncertainties in the observed 
Balmer edge depth play a major role in the uncertainty of 
inclination. For example, the two polarized spectra at wave- 
lengths shorter than about 300 nm are obviously below the 
model predictions. The spectral features were noticed by 
Kishimoto et al. (2004), who did not make unambiguous 
identification and interpretation of these features, but sug- 
gested that they might be related to Fe II absorption and 
possibly the Bowen resonance-florescence lines. Neverthe- 
less, the relatively lower signal to noise ratios of the mea- 
surements of the polarized spectra at wavelengths shorter 



than about 300 nm do not seem to affect the fitting results 
to the Balmer edge significantly. However, future better ob- 
servations may allow identification and physical modeling of 
these features, which should allow better determination of 
the Balmer edge. 



The foreseeable improvement in tightening the inclina- 
tion uncertainty is to have better polarized spectral observa- 
tions around the Balmer edge in the future. In this study, we 
find that the Balmer edge optical depth and the disc inclina- 
tion are particularly sensitive to the continuum emission at 
300 - 400 nm wavelengths. The current data have low signal- 
to-noise ratio in this wavelength range; with high signal-to- 
noise ratio data in the future, it is plausible to study the 
visible continuum of the quasar accretion discs, and con- 
strain the Balmer edge absorption and the disc inclination. 
Furthermore, since the uncertainty in inclination angle will 
propagate into the uncertainty in accretion rate for a given 
observed flux, a more accurate inclination will also result in 
a more accurate accretion rate. Therefore a combination of a 
more accurate inclination and a well observed peak emission 
will then allow the BH mass and accretion rate determined 
accurately, leading eventually to accurate BH spin measure- 
ment. 
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Figure 6. Correlation of the five fitting parameters of 4C 09.72. Each panel shows the joint posterior probability distribution between 
a pair of parameters with each dot corresponding to a parameter set in the MCMC chain of the second half of the 2 million iterations. 



APPENDIX A: ANALYTICAL EXPRESSIONS 

OF Kin, v AND R R 

We provide analytical expressions of the the (dimensionless) 
inner most circular orbit (x- ln ), the radiative efficiency (??) 
and the GR correction term (Rr), as functions of the BH 
spin (a/M) (Krolik 1998). The inner most circular orbit (xi n ) 
is 



X in = 3 + z 2 - x /(3-z 1 )(3 + z 1 + 2z 2 ) , (Al) 

for which z\ and z 2 can be computed from the BH spin 
(a/M) as 

z x = l+[l-(a/M) 2 ] 1/3 [(l+a/M) 1/3 + (l-a/M) 1/3 ] , (A2) 

(A3) 



~2 



yJz{a/MY + zl 



The radiative efficiency is related to the inner most circular 
orbit as 



xf n - 2xin + (a/M)x ; 



1/2 



(A4) 



The GR correction term (Rr) is a function of radius x 



Rr(x) 



C(x) 
A(x) 



(A5) 



where A(x) is 



A(x) = l- 2 - + ^l 

X X 2 



(A6) 



and C(x) is 
C(x) 



x yin 3(a/M) f „ 



v 2y 

3(2/i - a/M) 2 



yyi{yi - 


2/2) (2/1 - 


2/3) 


3(2/2 


- a/M) 2 




2/2/2(2/2 - 


2/i) (2/2 - 


2/3) 


3(2/3 


- a/M) 2 




2/2/3(2/3 - 


2/2X2/3 - 


2/1) 


VX, Vin 




and 



.'/in 

log- 
log 
log 



v - 


2/1 




- 2/1 


y - 


2/2 


Vin - 


-2/2 


y - 


2/3 


yin - 


- 2/3 



in which y 
roots of the following equation: 



x 3 -3x + 2a/M = 



(A7) 



(A8) 
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